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We investigate the electronic structure of SrsNiRhOs, a quasi-one-dimensional spin chain system 
using ab initio band structure calculations. Spin polarized calculations within GGA reveal that 
Ni and Rh have finite moments and they are antiferromagnetically coupled along the chain axis 
in the ground state. While these results obtained within the local spin density approximations 
provide remarkable representation of the magnetic phase, the experimentally observed insulating 
behavior could not be captured within this method. GGA+U calculations show that opening up of 
an insulating gap requires on-site Coulomb interaction among the Rh 4d electrons, « 2.5 eV and 
the correlation among Ni 3d electrons, C/J^' « 4.5 eV suggesting this system to be a Mott insulator. 
Electron correlation among d electrons leads to significant enhancement of the O 2p character in the 
energy bands in the vicinity of the Fermi level and the d bands appear at lower energies. Energy 
gap in the up spin density of states appears to be significantly small (~ 0.12 eV) while it is > 
2 eV in the down spin density of states suggesting possibility of spin polarized conduction in the 
semiconducting phase. 

PACS numbers: 71.27.-|-a, 71.20.-b, 75.20.Hr 



I. INTRODUCTION 

Recently, quasi-one-dimensional compounds, 
AEsMM'Oq {AE = alkaline earth metals such as 
Ca, Sr etc.; M and M' are transition metals) possessing 
rhombohedral K4CdCl6 structure (space group R3C) 
have attracted a great deal of attention due to their fasci- 
nating magnetic properties i^i^iii^i^iii^i^iiSiiiii^ii^iiiii^ii^ 
In this structure (see Fig. 1), MOe forms in trigonal 
prismatic geometry and M'Oq forms in octahedral 
geometry. These two units appear alternatively along 
the chain direction, c-axis and are connected via face 
sharing as shown in the figure. It is evident that the 
interchain interaction is significantly weak making 
these compounds a quasi-one-dimensional system. It 
is observed that antiferromagnetic, ferromagnetic or 
ferrimagnetic long-range ordered phases can be achieved 
by tuning the intrachain magnetic coupling in these one 
dimension chains i^i^ii In addition to such interesting 
properties of the quasi-one-dimensional chains, the whole 
system can be viewed as made off antiferromagnetic or 
ferromagnetic chains arranged in a triangular lattice. 
Numerous studies based on tailoring the composition of 
these geometrically frustrated systems reveal plethora 
of novel phases such as spin liquid phase, partially 
disordered antiferromagnetic phase eici^'^'^'^'^ 

In this paper, we report our results on the electronic 
structure of SrsNiRhOg using full potential linearized 
augmented plane wave method. This compound exhibits 
ferrimagnetic intrachain ordering below 45 K. Further 
lowering in temperature leads to partially disordered an- 
tiferromagnetic phase below 10 KA Interestingly, analo- 
gous compound SrsNiPtOg exhibits spin-liquid behavior. 
While it is difficult to capture disorder effect and/or spin 
liquid phase using such ab initio band structure calcula- 



tions, SrsNiRhOg is a good starting point to understand 
the intrachain coupling and associated ground state prop- 
erties. 

Our results clearly establish that the ground state of 
this compound consists of magnetic Ni and Rh ions which 
are antiferromagnetically coupled along the c-axis. To 
capture the insulating transport consistent with exper- 
imental observations, one needs to consider the on-site 
Coulomb interaction among Ni 3d and Rh 4d electrons 
in addition to the intrachain antiferromagnetic (lAFM) 
coupling, which suggests that this compound is a Mott in- 
sulator. On-site Coulomb interaction among Rh Ad elec- 
trons is found to be lower (about 60%) than that found 
for Ni 3d electrons as expected due to larger radial ex- 
tension of the Rh 4d orbitals. 



II. COMPUTATIONAL DETAILS 

The spin polarized GGA (generalized gradient approx- 
imation) and GGA-I-J7 {U = on-site Coulomb repulsion 
strength) electronic structure calculations were carried 
out using state-of-the art full potential linearized aug- 
mented plane wave (FPLAPW) methodJi The lattice 
parameters and atomic positions used in the calculations 
are taken from literature^. The Muffin-Tin sphere radii 
were chosen to be 2.33, 2.19, 2.01, and 1.78 a.u. for Sr, 
Ni, Rh, and O atoms, respectively. For the exchange cor- 
relation functional, we have adopted recently developed 
GGA form of Wu et al}^ The on-site Coulomb interac- 
tions were considered within LSDA-fJ7 (LSDA = local 
spin density approximation) formulation of Anisimov et 
al.^^ The calculations were performed by varying on-site 
Coulomb interaction among Ni 3d electrons, C/j^' and Rh 
4d electrons, U^j^. The convergence was achieved by con- 
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sidering 512 k points within the first Brillouin zone and 
the error bar for the energy convergence was set to be 
smaller than 10~^ Ryd/cell. 



FIG. 2: Color online) The partial density of states (PDOS) 
corresponding to do (thin solid lines), d±i (dashed lines) and 
d±2 (thick solid lines) orbitals for (a) Ni and (b) Rh obtained 
from spin-polarized GGA calculation when Ni and Rh are 
ferromagnetically coupled, (c) O 2p PDOS. 



III. RESULTS AND DISCUSSIONS 

The crystal structure of SraNiRhOg is hexagonal as 
shown in Fig. 1, where the Sr, Ni, Rh, and O ions are 
denoted by spheres of decreasing size. It is evident from 
the figure that the Ni and Rh ions are surrounded by 
six oxygen ions forming trigonal prism and octahedron, 
respectively. The trigonal prisms and octahedra are con- 
nected by face sharing and form the chains along the 
c-axis. If the axis system is defined such that z-axis is 
along c-axis and x- and y-axis are in the a6-plane, both 
Ni 3d and Rh 4(i orbitals defined in this axis system will 
not be degenerate due to the corresponding crystal field 
effect. There will be three bands consisting of do, d±i 
and d±2 orbitals. The energy separation and the loca- 
tion of the bands in the energy axis will depend on the 
type of crystal field is applicable as discussed later in the 
text. 

In Fig. 2, we show the partial density of states (PDOS) 



corresponding to Ni 3c?, Rh Ad and O 2p electronic states 
obtained from GGA calculations, where Ni and Rh are 
ferromagnetically coupled. The band splitting due to the 
crystal field as discussed above is visible clearly in the fig- 
ure. In Fig. 2(a), the energy range above -2 eV is essen- 
tially contributed by Ni 3d electronic states. The lower 
energy region (-4 to -2 eV) contains very low intensity of 
Ni 3d PDOS and has dominant O 2p contributions. The 
energy distributions of the density of states (DOS) corre- 
sponding to both Ni 3d and O 2p states appear similar. 
These results suggest that the DOS in -4 to -2 eV range 
can be attributed to bonding states arising from Ni 3c? 
and O 2p hybridizations and the antibonding states with 
dominant Ni 3c? character appear above -2 eV. 

The DOS in the energy ranges -6.3 to -4.4 eV and -1.5 
to eV correspond to Rh Ad-0 2p bonding and antibond- 
ing bands, respectively. The intensity of Rh 4c? PDOS is 
almost similar in these two regions, with negligible con- 
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tribution in the energy region of -4.4 to -1.5 eV (see Fig. 
2(b)). Two distinct features are observed in Fig. 2. (i) 
The energy separation between bonding and antibonding 
bands arising due to Ni 3d-0 2p hybridizations is signif- 
icantly smaller than the energy separation between the 
bonding and antibonding Rh Ad-O 2p bands, (ii) The 
bonding and antibonding Ni 3d-0 2p bands have domi- 
nant O 2p and Ni 3d character, respectively. However, 
the mixing of Rh Ad and O 2p characters in the Rh Ad-0 
2p hybridization is much stronger. These observations 
clearly demonstrate the stronger O 2p-Rh 4c? hybridiza- 
tion presumably due to larger radial extension of the 4d 
orbitals compared to the 3d orbitals. Such large over- 
lap integral, t leads to larger separation of the bonding 
and antibonding bands. However, the bandwidth of the 
individual bands does not appear to increase. 

In addition, it is evident from Fig. 2(a) that up-spin 
channel of Ni 3d states is almost fully occupied, whereas 
down-spin channel is partially occupied. Major contribu- 
tion in the unoccupied level comes from d±i orbitals. In 
the case of Rh Ad states, DOS corresponding to both the 
spin states are partially occupied and thus, contribute to 
the unoccupied DOS. 

The GGA calculations for ferromagnetic (FM) intra- 
chain coupling converged to the metallic ground state in 
sharp contrast to the experimentally observed insulating 
behavior in this compound. Moreover, the total magnetic 
moment per formula unit (fu) obtained from this calcu- 
lation is about 2.85 ^b, which is much larger than the 
experimentally estimated value of ^ 1 fiB-^ The mag- 
netic moment centered at Ni and Rh sites are 1.55 /is 
and 0.43 /is, respectively. Magnetic moment centered at 
the oxygen sites is found to be 0.13 /is, which is large 
and induced by the Ni 3d and Rh 4d moments. 

It is evident from the magnetic moments that if Ni and 
Rh are antiferromagnctically coupled, the total magnetic 
moment would match with the experimental results. In 
order to investigate such possibility, we have calculated 
the ground state energies and wavefunctions correspond- 
ing to antiferromagnetic coupling between Ni 3d and Rh 
4d moments. Interestingly, the calculated ground state 
energy is found to be 27 meV/fu lower than that of FM 
state. The calculated total magnetic moment comes out 
to be 1.04 /is/fu. Both these findings provide remarkable 
representation of the experimental magnetization datai^ 
The calculated magnetic moments for Ni and Rh are 1.31 
fj,B and -0.36 ^Bi respectively. The moment centered at 
oxygen site is almost negligible (0.02 /i^). 

The calculated Ni 3d, Rh 4d and O 2p PDOS are shown 
in Fig. 3. It is evident from the figure that intrachain an- 
tiferromagnetic interaction leads to significant reduction 
in bandwidth of all the energy bands and redistribution 
in spectral weight. For example, Ni 3d up spin band is 
almost completely filled and the Fermi level, is pinned 
at the top of the down spin band having dominant d±2 
contributions and a width of about 0.5 eV while Ni 3d 
down spin band in Fig. 2 has a width of about 1 eV and 
is partially filled. The occupancy of Rh 4d up spin and 
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FIG. 3: (Color online) The partial density of states (PDOS) 
corresponding to do (thin solid lines), d±i (dashed lines) and 
d±2 (thick solid lines) orbitals for (a) Ni and (b) Rh obtained 
from spin-polarized GGA calculation when intrachain Ni and 
Rh are antiferromagnetically coupled, (c) O 2p PDOS. 



down spin PDOS exhibit spectral distribution opposite to 
that observed in Fig. 2 as expected due to antiferromag- 
netic coupling. In this case also the bands are narrower 
than those observed in Fig. 2. The O 2p contribution 
at the Fermi level is dominant for up spin while in the 
ferromagnetic case it was down spin. 

All the above results clearly establish that the LSDA 
approach is adequate to capture the magnetic phase in 
this system. Although, there is significant spectral weight 
redistributions observed for different magnetic configura- 
tions, both these calculations converge to metallic ground 
state in contrast to the experimental studies. 

It is well known that such ab initio calculations often 
underestimate the on-site Coulomb correlation energy, 
which plays significant role in determining the electronic 
properties of d and / electron systems;^ The considera- 
tion of on-site Coulomb interaction among the d electrons 
(Udd) under GGA+U formulation is expected to improve 
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the situation. Thus, we have carried out GGA+U cal- 
culations for different values of Udd corresponding to Ni 
3d electrons, which is defined as /7j^* and Rh 4d elec- 
trons, which is defined as U^J^. The calculations were 
carried out considering the intrachain antiferromagnetic 
coupling. 

The calculated band gaps in the vicinity of for up 
spin and down spin density of states are shown in Fig. 
4. It is evident from the figure that consideration of the 
electron correlation for both Ni 3d and Rh Ad electrons 
helps to achieve the insulating phase. We have shown 
two types of results in the figure, (i) the dependence of 
the band gap as a function of U^J' for a fixed value of 
^dd" = 3 eV and (ii) band gap as a function of U^J^ for 
a fixed value of J7j^* = 5.5 eV. In the up spin channel 
(see Fig. 4(a)), it is evident that an increase in U^j^ does 
not lead to insulating phase till U^^^ = 4 eV. An increase 
in C/j^' above 4 eV provides an energy gap of about 0.12 
eV. Interestingly, the band gap remains unchanged for 
further increase in U^^\ On the other hand, the calcu- 
lations as a function of U^J^ for t/^' = 5.5 eV, indicate 
gradual increase in band gap. Thus, the energy gap is 
essentially determined by U^J^ when U^J" is kept large. 
The Fig. 4(a) suggests that the creation of insulating 
ground state requires U^J^ to be at least 2.6 eV. Since 
the Coulomb repulsion strength depends inversely on the 
separation of the two electrons, the strength of on-site 
Coulomb interaction is sensitive to the radial extension 
of the wave functions. The present result of [/^* larger 
than U^J^ is consistent with this behavior. 

Calculations for the down spin channel exhibit gradual 
increase in band gap with the increase in Udd for both Ni 
3d and Rh 4c? electrons. The magnitude of the gap in 
this case is always much higher (> 2 eV) than that (< 1 
eV) observed for up spin density of states. It is thus clear 
that various electronic properties of this compound will 
essentially be controlled by electronic density of states in 
the up spin channel. 

In order to investigate the character of the electronic 
states in the valence band, we show the calculated PDOS 
for U^/ = 5 eV and U^J" = 3 eV in Fig. 5. All 
the bands are significantly modified due to the spectral 
weight transfer from e^? to the energy region away from it. 
Ni 3d±i contributions exhibit significant band narrowing. 
The up spin bands appear around -5.8 eV and the down 
spin ones around 2.5 eV. The contributions from other Ni 
3d orbitals appear between -2 eV to -5 eV. Interestingly, 
the consideration of U^d leads to a significant increase of 
Ni 3d character at higher energies. Thus the O 2p con- 
tributions relative to Ni 3d become significantly high in 
the vicinity of the Fermi level. 

Consideration of the on-site Coulomb interaction 
among Rh 4d electrons reveals large effect on the density 
of states distribution. In the GGA calculation do was 
almost fully occupied and there was a significant con- 
tribution of d±i and d±2 orbitals within 0.28 eV of ep 
in the unoccupied up spin channel as evident from Fig. 
3(b). This trend is reversed in GGA+U calculations. 
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FIG. 4: (Color online) On-site Coulomb interaction among 
Ni 3d {U^J') and Rh 4d (U^Jl') electrons vs. insulating gap 
for (a) up-spin and (b) down-spin channels. Two cases are 
considered: (i) U^^ is fixed at 3 eV and J/J^' varied from 3.5 
to 5.5 eV (solid circles) and (ii) U^/ is fixed at 5.5 eV and 
U^^ varied from 3 to 5 eV (open circles). 

where region within 0.75 eV above the ep is contributed 
by the electronic states having do character. The elec- 
tronic states below ep are contributed by d±i and d±2 
orbitals. 

Despite the fact that the Ni 3d PDOS appear far away 
from ep, it has finite contribution in the vicinity of ep. 
This is evident in the rescaled Ni 3d PDOS (x 15 + 1) 
shown in Fig. 5(a) by dot-dashed line. Interestingly, the 
energy distribution of these density of states closely fol- 
lows the distribution observed in Rh 4d and O 2p PDOS. 
This clearly manifests the effect of hybridization between 
Ni 3d and Rh 4d states via O 2p states. Most interest- 
ingly, the gap in the up spin channel is much smaller 
than that found in down spin channel. Thus, this system 
provides an unique example of a semiconductor, where 
the electronic conduction is spin polarized. While half- 
metallic materials have long been investigated due to spin 
polarized conduction of charge carriers in this system. 
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FIG. 5: (Color online) The partial density of states (PDOS) 
corresponding to do (thin solid lines), d±i (dashed lines) and 
d±2 (thick solid lines) orbitals for (a) Ni and (b) Rh obtained 
from spin-polarized GGA+U calculation {U^J' = 5 eV and 
Udtt = 3 fiV) when intrachain Ni and Rh are antiferromag- 
netically coupled, (c) O 2p PDOS. 

This compound manifests properties that can be more 
useful in the semiconductor industry, where spin based 
technology is envisaged. 

Now, we discuss the manifestation of intrachain anti- 
ferromagnetic coupHng and Coulomb correlation in the 
band structure along different high symmetry directions 
of the Brillouin zone. The dispersion curves correspond- 
ing to different bands for FM, lAFM and lAFM con- 
sidering J7j^* = 5 and U^J^ = 3 eV for both the spin 
channels are plotted in Fig. 6. The zero in the energy 
scale indicates ep- In the ferromagnetic solution, the up 
spin bands along TM and TA directions do not cross the 
Fermi level, while the down spin bands cross ef in all the 
directions indicating half metallicity along these direc- 
tions. The antiferromagnetic coupling among Ni and Rh 
in the chains leads to narrowing of the bands and the up 
spin bands cross the Fermi level in all the directions. On 
the other hand the down spin bands in the vicinity of ep 



FIG. 6: Dispersion curves along different symmetry directions 
of the first Brillouin zone for both the spin channels. Fermi 
level is denoted by zero in the energy scale. 



essentially appear below ep. Again, the band crossing of 
the down spin band is observed only along AT direction 
suggesting half metallicity in other directions. 

Inclusion of the electron correlation among the d elec- 
trons leads to significant change in the band structure. 
While, there is a hard gap in both the spin channels con- 
sistent with the transport behavior of this compound, 
the gap exhibit significant fc-dependence. The gap at the 
high symmetry point, F is lower than that observed at M, 
K or A points. Interestingly, the dispersion of the bands 
close to ep appears to be larger than that observed in 
the uncorrelated case. This indicates that although elec- 
tron correlation essentially tends to localize the electrons, 
the electronic states in the vicinity of the Fermi level are 
highly mobile. 

The occupancies of different d orbitals corresponding 
to Ni and Rh obtained from GGA+U {U^J = 5 and U^J" 
= 3 eV) calculation are shown in Table 1. In up spin 
channel, all the Ni 3d orbitals are almost completely filled 
and in down spin channel, only do orbital is completely 
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filled. The small deviations (~ 0.05) from the completely 
filled values can be attributed to the hybridization of Ni 
3d orbitals with the neighboring O 2p orbitals. Total 
number of electron in the Ni 3c? is found to be 7.91, which 
suggests that the valency of Ni is close to 2+. 



Table 1: The electron occupancies of , d±i, and 
d±2 orbitals for Ni and Rh obtained from spin polarized 
GGA+J7 calculation. 





d-2 


d_i 


da 


di 


d2 


Total 
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0.95 


0.98 
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0.98 
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Ni-down 


0.8 


0.28 
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3.11 


Rh-up 


0.65 


0.5 


0.21 


0.5 


0.65 


2.51 


Rh-down 


0.66 


0.52 


0.82 


0.52 


0.66 


3.18 



In the case of Rh 4d orbitals, both the spin channels 
are partially occupied. Our calculations exhibit 5.69 elec- 
trons in the Rh 4c? orbitals. If we consider O in 2- and Sr 
in 2+ states, then charge neutrality condition demands 
Rh to be in 4+ state, which corresponds to 5 electrons 
in the 4c? orbitals. The extra 0.69 electrons can be at- 
tributed to strong hybridization of Rh Ad orbitals with 
the neighboring O 2p orbitals. 

The calculated magnetic moments for Ni and Rh ions 
are 1.69 fiB and -0.68 fiB^ respectively, with almost neg- 
ligible contribution from oxygen (0.01 /is). The total 
magnetic moments per formula unit, including the contri- 
bution of -0.07 iiB from the interstitial region, is found to 
be about 1 /is consistent with the experimental results.^ 

In order to understand the possible cause for lAFM 
interaction between Ni and Rh spins, we show the 
schematic crystal field split energy level diagram for both 
Ni and Rh obtained from spin unpolarized GGA calcula- 
tion by considering centers of gravity of c?o, c?±i, and c?±2 
PDOS is shown in Fig. 7. It is evident from the figure 
that c?o orbital has lowest energy for both Ni and Rh. 
However, next energy levels consist of different orbitals. 
In Ni, next level corresponds to d±i orbitals, whereas 
in Rh it is c?-|-2 orbitals. The separation of Rh d±2 and 
c?±i is found to be ~0.7 eV, which is about double of 
the Hund's couphng energy (Jh) of Rh 4c? electrons es- 
timated from spin polarized GGA calculation. On the 
other hand the energy separation of Ni c?o, c?±i, and d±2 
orbitals are found to be ^0.3 eV, which is about one fifth 
of the estimated value of Jh for Ni 3c? electrons. 

As discussed above, Ni is in 2+ state leading to 5 = 
1 state. Rh'^+ has 5 electrons in the 4c? orbitals. The 
schematic diagram clearly shows that Rh atoms will be 
in low-spin (5=1/2) state with 2 electrons occupying c?o 
orbital and the other 3 electrons will occupy c?±2 orbitals 
making them more than half-filled. It is well known that 
in such a case, the interacting electrons via oxygen 2p or- 
bitals (superexchange interaction) favor antiparallel cou- 
pling in the ground statej^ 



dj d 




FIG. 7; (Color online) Schematic energy level diagram of do 
, d±i, and d±2 orbitals for Ni^"*" and Rh*+ showing antiferro- 
magnetic coupling as the possible ground state. 



IV. CONCLUSIONS 

In conclusion, we studied the electronic and magnetic 
properties of SraNiRhOg using ab initio spin polarized 
band structure calculations within GGA and GGA+U. 
The calculations for intrachain antiferromagnetic cou- 
pling exhibit the lowest energy. The magnetic moments 
of Ni and Rh ions were found to be '^1.7 /is and -0.7 
fiB, respectively, with a total moment of about 1 fiB- 
All these results corresponding to the magnetic phase 
provide a remarkable representation of the experimental 
results. 

The insulating phase observed in transport measure- 
ments could not be captured within GGA suggesting im- 
portance of electron correlation. GGA+U calculation 
reveals that the critical values of Coulomb correlation 
strength to open up the insulating gap are ~ 4.5 and 
~ 2.5 eV for Ni 3d and Rh 4d electrons, respectively. 
The inclusion of electron correlation leads to significant 
change in character of the energy bands. The feature 
close to ep exhibits enhancement of O 2p character, while 
the d character enhances at higher energies. Interestingly, 
the band gap in up spin channel is found to be signifi- 
cantly small (^ 0.12 eV) compared to that observed (> 
2 eV) in the down spin channel indicating possibilities of 
spin polarized conduction in the semiconducting phase. 
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